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1. Introduction 

In the past few years we have learned how to generate superconformal theories in 
various dimensions. Nevertheless we do not understand the properties of conformal theo- 
ries in dimensions greater than two with the same precision as we do in the case of two 
dimensions. 

In two dimensions, one powerful method to study conformal theories was pioneered 
by Zamolodchikov [jl|, who studied their properties under special integrable massive de- 
formations. Zamolodchikov's idea was to use the properties of the deformed theories, such 
as the spectrum of kinks and their scattering matrix amplitudes, to reconstruct the full 
data of the conformal theory. As a particular example of this philosophy, it was subse- 
quently shown in the case of massive deformations of (2, 2) superconformal theories in two 
dimensions 0] that just the degeneracy of the (BPS saturated) kinks of the massive theory 
characterizes the dimensions of chiral operators at the conformal point. 

It is thus natural, in the context of conformal theories in dimensions bigger than 

2, to also ask the same question, namely, what is the number of BPS states for slight 
deformations of these theories away from the conformal point. For certain cases, this has 
been answered. For example we know the degeneracy for A^ = 4 theories in 4 dimensions 
and some superconformal N = 2 theories (for example SU{2) with four doublets) [Q. The 
aim of this paper is to widen the class of theories for which we know the BPS spectrum. 

The case of A = 2 field theories in 4 dimensions is already very interesting from this 
point of view. In a limited set of cases, including certain A = 2 SU(2) gauge theories, the 
BPS spectrum can be computed exactly using field theory techniques 0] 0. However, there 
are many interesting 4D A = 2 theories for which we know very little about the structure 
of their BPS spectra, including the superconformal theories of Argyres and Douglas 0. In 
this paper, our main result will be a determination of the degeneracy of the BPS states un- 
der slight deformations away from Argyres-Douglas points. Q We realize Argyres-Douglas 
superconformal theories and their deformations in terms of type IIB strings propagating 
in a (nearly) singular Calabi-Yau threefold. The BPS states are realized in this setup in 
terms of D3-branes wrapped around supersymmetric 3-cycles on the threefold. The struc- 
ture of BPS states of Argyres-Douglas superconformal theories are among the simplest to 
study from this point of view because the counting of the supersymmetric 3-cycles gets 



A previous study of the spectrum of an Argyres-Douglas point using M-theory has been 
made by Gustavsson and Henningson 0. 
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related to counting certain special 1-cycles on a Riemann surface. We find a structure very 
reminiscent of BPS degeneracies of the deformations of minimal N = 2 superconformal 
theories in 2D. In particular we find that the spectrum of light BPS states is finite and 
depends on how we deform away from the conformal point. 

The structure of this paper is as follows: In section 2 we review certain aspects of 
Calabi-Yau (i-folds which have isolated singularities. We point out a connection between 
normalizable deformations of a Calabi-Yau singularity and unitarity bounds of the corre- 
sponding quantum field theory in the remaining R"^. In section 3 we make some remarks 
about the string worldsheet description near a Calabi-Yau singularity. In section 4 we 
specialize to the case of Calabi-Yau threefold, which gives rise to A^ = 2 superconformal 
theories in 4 dimensions. We set up the computation for counting the number of super- 
symmetric 3-cycles for this case. In section 5 we apply these techniques to the case of 
Argyres-Douglas conformal theories and compute their spectra. 

2. Calabi-Yau Singularities 

In this section we consider Calabi-Yau manifolds with an isolated singularity. More 
precisely we consider a local model for a d-iold given by a hypersurface W{xi) = in C*^"*"^. 
Moreover we assume W{xi) is a quasihomogeneous function of Xii 

W{X'i'Xi) = XW{x,). (2.1) 

We assume the singularity is isolated, namely 

dW = iff Xi = 0. 

If W is viewed as a LG superpotential of an A^ = 2 theory in 2 dimensions it would 
flow to a superconformal theory with central charge c given by 

d+l 

c = 5^(1 - 2q,) = {d+l)- 2(^ g,)- (2-2) 

i=l i 

The condition that the Calabi-Yau threefold with this singularity appear at finite distance 
in moduli space was recently studied in |T^ with the conclusion that if 



c<d-l (2.3) 
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the distance is finite and thus the study of singularity would be of physical interest only 
in this case. 

An interesting question involves deformations of the singularity. In this context, as in 
the 2 dimensional case [§] , it is natural to consider the singularity ring of W 

7^ = C[xi]/dW 

generated by the monomials 

^Oi _ ai 02 i^d+l 

•Xj tAj -t %Li ty • • • Uj J_|_ 1 

modulo setting to zero those polynomials which are in the ideal generated by di^ . One 
then considers deformations of the form 

For each element x" in the ring, consider the charge 

i 

The Qot lie in the range 

< Qa < c. 

Moreover, the elements in the ring are paired so that for every element of the ring with 
charge Qa there is one with charge c — Qa- The degeneracy of the ring elements is captured 
by the function 

and the dimension of the ring Tt is given by 

d+l 



A^ = dim7^= JJ 



(1 - ^^) 



i=\ 



This dimension is also the dimension of compact part of mid-dimension homology H^iy^ = 
\i) . The compact homology can be realized by N spheres of dimension d with an intersection 



structure that can be obtained from the structure of VK [pT . 

In questions in string theory, it is important to know how many cohomology elements 
of dimension d are supported on the singularity. These would correspond to normalizable 



d-forms in the internal dimensions localized near the singularity. These forms can be 
obtained from the holomorphic d-form represented as 



dW 



What this means is that we solve the equation 1^ = for one of the Xi in terms of the 
others and view dxi/dW as 1/diW in the above expression. The other d-iorms correspond 
to deformations of O: 

Qa = dfl/dga- 

The condition that Q^ lead to a normalizable form localized at the singularity was studied 
IPf . The condition for this is that J |Oq,P diverge in the vicinity of Xi = 0. It was found 
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that the deformation by monomial Xa corresponds to a normalizable cohomology element 

Q. < ^^. (2.4) 

Note that combined with ( |2.3D this implies that a normalizable cohomology deformation 
must have Qq < 1. 

The condition ( p.4|) can also be understood physically by unitarity bounds of physical 
modes as follows. The coefficients g^ in the deformation of the singularity should corre- 
spond in the field theory setup, either to expectation values of scalar fields, or to coupling 
constants of field theory. We may assign dimensions to the coefficients g^ as follows: If 
we consider a Dd-hrane (or a d-dimensional part of a higher dimensional brane) wrapped 
around a (i-dimensional supersymmetric cycle C in the local geometry, it gives a particle 
with BPS mass 

Jc 
which implies that the mass dimension of [O] = 1. By quasihomogeneity, we can assign 
a dimension to each variable Xi proportional to its charge. The condition that [O] = 1 
implies that the dimension of a monomial x^ is given by 

and that the dimension of g^ is given by 

r 1 2(1 -Q^) 

[9.] = j^zrr^y (2-5) 
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The deformation parameter Qa is the expectation value of a canonicaUy normahzed scalar 
field weighted with the Yang-Mills coupling; that is, Qa = gYM{4>a)- Since [(/^m] = 
(4 — -D)/2 in D spacetime dimensions, the unitarity bound that [(j)a\ > (-D — 2)/2 in a 
superconformal theory requires that 

r 1 2(1-0^ 



and this condition is identical to ( p.4| ). Thus, the mode is a physical field if and only if 
it satisfies the unitarity bound. This is quite a satisfactory result (which has also been 
obtained independently in [|T2|). 

Let us now consider Calabi-Yau manifolds of various dimensions and see what this 
condition translates into in each case. 

2.1. d=2 

In the case of complex dimension 2, the finite distance condition (J2.3D c < 1 implies 



that W is given by the LG superpotential of A^ = 2 minimal models, which gives the usual 
ADE classification of K3 singularities. Moreover the condition of physical fields ( p.4| ) is 
automatically satisfied for all elements in the ring because all the elements have 

Qa<C< ——. 

This is in accordance (in the type IIA case) with the fact that all these deformations can 
be viewed as giving vev to scalars in commuting directions of the adjoint representation of 
the corresponding group. 

2.2. d=3 

For Calabi-Yau 3-fold singularities, the condition of being at finite distance implies 
that c < 2. If we consider type IIB strings in the presence of such singularities we obtain 
an A^ = 2 superconformal theory in 4 dimensions. Examples of such 1^'s are provided by 

W = F{x,y) + z'^ + w'^ (2.6) 

where F{x,y) is a quasihomogeneous function of x,y. In this case the function F{x,y) 
can be identified WM with a fivebrane with worldvolume R^ x "E, where E is the Seiberg- 



Witten Riemann surface F{x, y) = 0. In this case the singular Calabi-Yau gets mapped 
to a singular Riemann surface. A special case of this 



n I „,2 



F{x,y) = x'' + y 



corresponds to the Argyres-Douglas points that we will study in more detail later in this 
paper. 

There are other singularitities which are also at finite distance which are not of the 
form (|2.6| ). For example consider 

Type IIB string in the presence of such a singularity describes a superconformal theory, 
which does not have a description in terms of a Seiberg-Witten Riemann surface. Such sin- 
gularities do arise in certain gauge theories. For example the above singularity appears in a 
gauge theory with gauge group SU{3N) x SU{2N)'^ x SU{N)^ with certain bi-fundamental 
matter dictated by the affine Eq Dynkin diagram ||14|| . 

In the case of the threefold, the condition that deformations x" correspond to nor- 
malizable forms is ( |2.4| ) 

Qa < c/2. 

Given the pairing of the ring elements and the fact that the sum of the charges of the 
pairs gives c, this means that the deformations that correspond to expectation values of 
dynamical fields are in one to one correspondence with the deformations that correspond 
to parameters in the theory. A natural interpretation of this pairing is to note that for 
each chiral operator in the N = 2 superconformal theory $q, (whose lowest component is 
identified with (p^ where {(pa) = ga)^ we can consider the N = 2 superspace integral 

s ^ s+ f d^xd^eta^a 

with tct being identified with the dual deformations of W. The form of these deformations 
implies that the mass dimensions obey [ta] + [Qa] = 2. In fact this is consistent with what 
we have found for the dimension of the parameters. Namely, since the charges of dual 
deformations add up to c we deduce that the charge corresponding to the ta deformation is 
c—Qoi and using the dimension of the fields given by (|2.5|) we see that indeed [ta] + [fi'a] = 2: 

r 1 2(1 -gj 



[tc 



(2-c) 
2(l-c + Q„: 



(2-c) 
[ga] + [U = 2. 
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2.3. d=4 



The case of local singularities of Calabi-Yau 4-folds was studied in |T^. In particular 



the case of theories with c < 1 was studied in detail and it was shown that type IIA string 
in the presence of these singularities gives rise to certain N = 2 superconformal Kazama- 
Suzuki models in 2 dimensions. In these cases none of the deformations correspond to 
dynamical fields in the 2d theory. 

3. Perturbative String Description 

It is natural to consider the perturbative string theory in the presence of such singu- 
larities. Of course, if we are at the singularity the string theory is singular and there is 
no perturbative expansion. However, if we deform the quasihomogeneous function W we 
can resolve the singularity. In such a case (if the singularity is resolved enough so that the 
wrapped branes are heavy enough) we can expect to have a perturbative string description. 
In this case one can use ideas developed in |]l^ [|l6| in the context of compact Calabi-Yau 



manifolds and for the non-compact case in [|T^ [^ |]T^ to find a Landau-Ginzburg descrip- 
tion of the theory. Namely consider d + 2 chiral fields xi, .., x^+i, y with N=2 U{1) charge 
given by gi, ..., Qd+i, —1/7, and with a quasihomogeneous LG superpotential 

Vr = Vr(xi)+^^ax"y^('5"-^). (3.1) 

ct 

7 is fixed by the requirement that the total c = d, which implies that 

2 
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d-c-1 

Note that the condition that the singularity be at finite distance ( |2.3| ) is that 7 > and 
the condition that the deformation given by g^ corresponds to a normalizable field in 
uncompactified theory is translated to the condition that the power of y accompanying 
that deformation be less than —1. To obtain the full supersymmetric theory we also need 
to do the Gepner projection, which is a discrete orbifold of the above LG theory that keeps 
only the integral U{1) charged fields. 

A particular case of ( P?T| ) is the deformation of VF by a constant: 



In this case it can be shown [^ |T^ [jT^ that the resulting theory corresponds to the tensor 



product of the LG theory given by W and a Kazama-Suzuki coset model SL{2)/U{1) with 
level k = 2 -f 7. This construction has also been noted independently in [O. Certain 



aspects of conformal theories near singularities have also been discussed recently in pO 



4. Solitons in the d = 3 case 

For the rest of this paper we will be interested in studying some aspects of Calabi-Yau 
threefold singularities corresponding to Argyres-Douglas points. As already noted these 
points correspond to local singularities given by 

W = P{x) + y'^ + z'^ + w^ = (4.1) 

with P{x) a polynomial of degree n in x. The superconformal point itself corresponds to 
P{x) — x"', and there are n — 1 deformations given by 

n-2 

P{x) = x" + ^ gix' 

i=0 

where the gi correspond to expectation value of fields for i < (n — 2)/2 and to dual mass 



parameters for i > (n — 2)/2 [21|. Eq.( |2.2| ) gives the central charge 

2 



1- 

n 



and (^^) gives the dimensions 

. . 2(1 -z/n) 2(n 



1 + 2/71 n + 2 



0..n-2) 



in agreement with the result of [p2 



We are interested in studying type IIB strings in the presence of such singularities. In 
particular, we would like to study the BPS states in this theory near the superconformal 
point and how their spectrum jumps as we change the parameters gi in the polynomial 
P{x). 

The BPS states in this case correspond to D3-branes wrapped around supersymmetric 
3-cycles. A basis of vanishing 3-cycles (not necessarily supersymmetric) can be chosen to 
be n— 1 3-spheres intersecting one another according to the Dynkin diagram of A„_i. The 
intersection of two cycles can be interpreted as the skew symmetric form in the product of 
electric and magnetic charges. In particular, two S^^s corresponding to adjacent nodes of 
An-i carry, in some basis, electric versus magnetic charge. 



To find the supersymmetric 3-cycles, one follows the strategy in |T^ and considers the 
3-cycle as a 3-sphere consisting of an 5''^ fibered over a real curve in the x-plane. The 5''^ 
above a given x is y"^ + z"^ + w'^ = —P{x) (in an appropriate real subspace). The projection 



of the 3-cycle onto the x-plane is a curve in the a;-plane which begins and ends at zeroes of 
P{x). At the endpoints of the curve, which correspond to the 'poles' of 5"^, the radius of 
the S'^ in the fiber goes to zero. A basis of 3-cycles (not necessarily supersymmetric) can 
be chosen to correspond to n — 1 intervals connecting the n zeroes of P{x) in a sequence. 
These would correspond to a set of 3-cycles whose intersection gives the Dynkin diagram 

of An- 1. 

As discussed in [|15| the condition of having a supersymmetric cycle gets translated in 
the x-plane into the existence paths beginning and ending at the zeroes of P{x) such that 
the phase of 

[ n = ^/P[x) dx (4.2) 

is constant along the path. This would guarantee in particular that the BPS inequality 



|0| > I / 0| 
'S3 Js^ 

is saturated. Alternatively, the required condition is that the image of the path under the 

Jacobian map with respect to the reduced one-form should be a straight line in the flat 

VK-plane, that is, 

W{x{t)) = / VP(x) dx = at (4.3) 



where t is real parameterizing the path and a is a phase. In this case, the mass of a D3- 
brane wrapped around the supersymmetric cycle corresponding to such a path between 
zeros at xq and xi automatically saturates the BPS inequality 



M 



Xl 



Xo 



Pdx 



> 



Xl 



Pdx 



xo 



An efficient technique to find such paths is to solve the first-order differential equation 

dx 



^/P{x) 



dt 



a. 



for solutions x{t) beginning and ending at roots of P{x), for all possible values of the phase 
of a. This is equivalent to finding integral curves of the vector field 

" A (44) 

that start and end at roots of P. As we have explained, the existence or nonexistence of 
such a curve then implies the existence or nonexistence of a particular BPS state. This 
procedure is easily implemented. 

In the next section we analyze the solutions to the condition ( [4.3| ) in various regimes 
of parameters. 
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5. Finding the BPS states 

As we have discussed, we need to find solutions to 

dx a 



dt ^/P(x) 



(5.1) 



where t is a real parameter and a is a phase. Moreover the solution should begin and end 
at one of the n roots of P{x). We need to establish a few facts: 

For each topology of path between two roots of P{x), a is uniquely fixed (up to an 
overall sign). By topology of path we mean a homology class of paths with fixed endpoints 
on the complex x-plane with all the roots of P{x) deleted. To see the uniqueness of a note 
that for each path 7 in a particular homology class the integral 

dx^/p(x) = a^. (5.2) 

which is well defined up to an overall sign, is independent of the precise choice of 7 in 
that class. This is because two different paths in the same class will differ by integral of 
an analytic expression in the region bounded by the two curves. Physically, this is simply 
the refiection of the fact that each choice of path topology fixes the electric and magnetic 
charges of the BPS state and the integral ( |5.2| ) is just the central charge of the N = 2 
algebra in that sector. It follows that a is uniquely fixed for each choice of class of path 
to be 

a = a.j/\a.j\. 

The next fact we need to establish is that there is at most one solution connecting two 

roots. We first show that there cannot be two different solutions in the same path class. 

Let us first analyze the structure of the solution to ( [5.1| ) near each root. Near a root of 

P{x), which with no loss of generality we take to be at x = 0, one is solving an equation 

which can be approximated as 

dx a 

dt -Jx 



whose solutions are given by 



X — (|at) 



In particular, for a given a. there are three solutions near 3: = which make 120° angles 
relative to one another, corresponding to the 3 possible choices of cube roots of o? in this 
solution. A typical set of integral curves near to a root is depicted in Figure 1. 
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Fig.l; A typical set of integral curves near to a root. 

Now suppose there are two solutions in a given class (and thus with the same choice 
of a) connecting two roots of P. In other words, let us assume that two of the three 
integral curves from one root join with two of the three integral curves from another root 
(see Figure 2a). The geometry should be as depicted in Figure 2a: first, the paths cannot 
intersect each other except at the two ends, as the vector field has a well defined direction 
at each point and that would not be the case at an intersection point. Secondly, the paths 
must make an angle of 120° on each side, otherwise the third integral curve emanating 
from either of the roots would have nowhere to go. (We can also rule out the case where 
all three curves are joined, by applying this argument to two of the curves making 120° 
angles relative to each other.) From the geometry of Figure 2a it is clear that there will 
necessarily be closed integral curves trapped inside. If we consider the winding of the phase 



defined by 



$ = lv/pM|/v/pM 



w = / d log $ 

27rz 



we see that w = +1 along these curves. Since P{x) has no zeroes or poles inside this curve, 
this is impossible. We could also have chosen, instead of the trapped integral curve, an 
arbitrary interior curve to use for counting the winding of the phase $. For example we 
can make the choice shown in Figure 2b. In this case we again get winding number +1, 
with the two small arcs near the roots each contributing +1/6 to the winding number and 
the remaining pieces along the integral curves giving a net winding of +2/3. 
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Figs. 2a and 2b: Impossibility of more than one integral curve of a given topology between two 
roots. The winding of the phase is calculated by integrating along the dashed line. 



Now let us also establish that there can not be two solutions between two roots, even 
if the paths are in different topological classes. In such a case the two paths correspond 



to solutions to (|5.1|) with two different phases ol\ and ol2- Let e = ol\ol^ and define 
b = 2^ log e. With no loss of generality we can assume e to be some generic phase (in 
other words for any special values of e we can change the coefficients of P{x) so that the 
two central charges have infinitesimally different phase ratios). Let us consider the winding 
of the phase $(x) along the closed path shown in Figure 3. We will first assume, as in 
Figure 3, that the two curves do not intersect except at the endpoints. This assumption 
will be justified shortly. 




a, 



Fig. 3; Two integral curves between roots, with ai ^ ai- 
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This closed path consists of the two integral curves and the small paths connecting them 
near the roots. Even though the two integral curves are for different choices of a, in either 
case a is a constant and does not vary over the path. In particular the integral J d log $ 
along the two integral curves is the same as winding of the velocity vector, and gives a total 
contribution to the winding of +2/3. Along the arcs near each of the roots, the contribution 
to the winding will be given by 1/6 up to a correction by 6/6 at one endpoint and —6/6 
at the other endpoint, which cancel out. We thus still obtain a net winding of +1 unit. 
However this is impossible as the winding associated with 1/Vp is in general negative and 
given by —m/2 where m denotes the number of roots of P inside the closed curve. The 



reason for this is that 1/y/x — a has a winding number —1/2 around x = a. (The other 
possible geometries for these two integral curves curves — in which their relative angle of 
approach to one of the endpoints is shifted by ±2n/3 — lead to similar contradictions.) 

a 




«2 

Fig. 4; Impossibility of intersection of two integral curves. 

We also need to justify our assumption that the curves do not intersect each other 
except at the endpoints. Suppose they did. Pick an endpoint and consider the closed curve 
that forms between it and the first point where the two integral curves intersect (Figure 
4). The winding number around that closed curve is again the sum of contributions of the 
two integral curves and the small arcs near the sharp ends. The integral curves contribute 
the sum of the two angles divided by 27r, or ^(1 + 5) + 5, the left endpoint adds |(1 + 5), 
and the right endpoint adds nothing. So the net winding number | + |5 is not even an 
integer, as e can be assumed to be a generic phase. This is clearly a contradiction with 
the fact that the winding should be —m/2 for some integer m. 

We have thus established that there is at most one solution connecting two distinct 
roots. 

Minimum Number of BPS states: We can also argue that any pair of roots are con- 
nected by a sequence of BPS solutions. We can find such a sequence by minimizing 
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J \^^/P{x)dx\ over arbitrary paths beginning at one root and ending on the other. Clearly 
there will be a solution to the minimization condition (paths going to infinity will give an 
infinite contribution to the above integral). Consider the path which minimizes it. Then 
for a generic point on this path, away from the roots of P, the curve should satisfy ( |5.1| ). 
If not, fix two nearby points on the curve and take a solution of ( |5.1|) between these two 
points. This will have a lower value for the integral, violating the assumption that we had 
found the curve which minimizes the integral. In general, this minimal curve may pass 
through some roots of P. The above argument still shows that piecewise, between the 
roots, it should satisfy ( p.l| ) for some a. This implies that any two roots are connected 
by a sequence of BPS solutions. If there are n roots in all, the minimum number of BPS 
solutions is thus n—1. 
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(a) (b) 

Fig. 5a: Integral curve from 1 to 3 when 9 < 120°. Fig. 5b: When 9 > 120° no 1-3 curve exists. 



Jumping Phenomena: Consider three roots of P labeled by 1, 2, and 3. Suppose for 
some choice of parameters there are BPS solutions from 1 to 2, from 2 to 3, and from 1 to 
3. A priori, the phase of the central charge of the N = 2 algebra, which determines the 
phase a, will be different for the three different solitons. However, let us assume that as 
we change the parameters the phases become the same, i.e., the BPS charges from 1 to 
2 and 2 to 3 get aligned and sum up to the charge from 1 to 3. This means that for the 
same choice of a, i.e. ai2 = 023 = ciis, we have three integral solutions. In such a case the 
curve from 1 to 3 must coincide with the concatenation of the curves from 1 to 2 and 2 to 
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3. If not, as in the previous arguments, by considering the trapped integral curves inside 
we would get a contradiction with the winding number. At the alignment point the angle 



= 123 



which gives the angle between the solitons 1-2 and 2-3 is 120 degrees. However as we 
change parameters from one side to the other side the angle changes from less than to 
greater than 120°. We will argue below that the case where the 1-3 soliton does exist 
corresponds to when 6 < 120°). First, we argue that the 1-3 soliton should disappear 
when 9 > 120°. If the 1-3 solution continued to exist, the case before (when 9 < 120°) 
and after (9 > 120°) the transition should look like that depicted in Figures 5a and 5b. 
This is because the phase of Q;i2/ai3 should go from one sign to another, which implies, by 
looking at integral curves near the point 1, that the 1-3 curve is to one side or the other 
of the 1-2 curve. But the curve 1-3 depicted in Figure 5b is forbidden: The third integral 
curve emanating from point 2 would intersect it, which is not allowed. We thus see that 
the soliton represented by the 1-3 curve decays to 1-2 and 2-3 soliton as we pass through 
alignment. 




1 3 

Fig. 6; Proof of existence of integral curve from 1 to 3 when 9 < 120°. 

The reverse of this argument can also be made. Suppose we originally have no soliton 
from 1 to 3. Let us choose parameters such that we are close to an aligned configuration, 
which as argued should correspond to 6* = 120°. Then we argue that as we decrease 9 a 
soliton should appear from 1 to 3. To see this note, as we discussed already, that there 
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is a solution which minimizes the integral J \y/Pdx\ between any pairs of point. Apply 
this to points 1 and 3. At the alignment point the minimal solution is given by the sum 
of the curves 1-2 and 2-3. If 9 is decreased to less than 120° the class minimizing it is 
no longer sum of the 1-2 and 2-3 curves. The integral can be decreased, for example by 
considering a path such as the one shown in Figure 6: we pick two points A and B, on 1-2 
and 2-3 curves respectively, near root 2. We consider a solution of the integral curves •jab 
of the vector field which interpolates between A and B for some a. That this is possible is 
guaranteed if ^ < 120°. Then the integral J \y/P dx\ is decreased if for the part near point 
2, the curve •jab is used, instead of the radial lines connecting A and B to point 2. Thus 
the minimum does not coincide with the sum of the 1-2 and 2-3 curves, and there appears 
a distinct 1-3 soliton as in Figure 5. The local analysis here is similar to what Joyce has 
considered in the context of 3-cycles on Calabi-Yau [^ Q. 

To summarize, whenever the phases of 3 such BPS charges are equal, we are in a 
situation of marginal stability. Under a perturbation away from alignment in one direction, 
all 3 states will be stable, but in the other direction, one of the states will become unstable 
to decay into the other two states. 

5.1. Examples 

In this part we show that with suitable choices of P{x) we can have any number of 
solutions between n — 1 and n{n — l)/2. That ti — 1 is the minimum possible number has 
already been noted. Also that n{n — l)/2 is the maximum number follows from the fact 
that there is at most one BPS state for each pair of roots of P. We will show that, for 
example, with P{x) = x"^ — 1 we get n{n — l)/2 solitons, one for each pair of roots of P. 
On the other hand with P{x) having only real roots, we show that P has exactly n — 1 
solitons, corresponding to the solitons connecting adjacent roots. By our discussion about 
the jumping phenomenon, it follows that as we continuously change the polynomial P{x) 
we get arbitrary number of solitons anywhere between these two bounds. In this way the 
story is very similar to that of the An series for N = 2 LG theories in 2 dimensions eI. 

Let us first focus on the simplest nontrivial A^ singularity, with P{x) = x^, corre- 
sponding to the original SU(3) Argyres-Douglas point p. By the above arguments, the 



We would like to thank M. Douglas for pointing this reference out to us. See also [24| [25|. 
In fact if we consider a degenerate choice of polynomial, where P = (dW/dx)'^ the problem 



of solving (5.1) is identical to the problem of finding solutions to 2d solitons in an A'^ = 2 theory 



with LG potential given by W. 
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miniinuni possible number of BPS states is n — 1 = 2, and the maximum is 3. A general 
deformation of P(x) 

P{x) — X^ + UX + V 

involves two parameters u and v, with v an operator expectation value and u a mass 
parameter. As we shall describe, the stable BPS spectrum near the singularity depends 
on the direction of approach, and the marginal stability surfaces (MSS's) extend all the 
way down to the critical point. By tuning the dimensionless ratio xc' jv^ appropriately as 
the scaling is performed, we may end up with either two or three light stable BPS states. 
(This modifies the result of 0, which found three states.) We will now exhibit examples 
in which each of these possibilities is realized. 

First consider a situation where all three roots are arranged symmetrically, at cube 
roots of unity times a common scale factor, so that P{x) = x^ — 1. Given any integral 
curve linking two roots at some value of a, we can construct two others, linking the other 
two pairs of roots, by multiplying a. and x by an appropriate cube root of unity (which 
preserves P{x) and the form of Eq.( ^?3|) ). Hence, the total number of integral curves must 
be at least 3. Since the total number of such curves must also be either 2 or 3, this 
maximally symmetric configuration of roots leads to exactly 3 stable BPS states. 

As another example, suppose the roots are all colinear; without loss of generality, we 
may suppose that 

P{x)^x{x-\){x-\) (5.3) 

where A is real. We claim that there are exactly two integral curves for all A > 1. By 
symmetry, all integral curves must either be at a = ±1 or a = ±z, and must lie along 
the real axis. Otherwise, given an integral curve not satisfying these conditions, we could 
construct another inequivalent geodesic between the same two roots, at the conjugate value 
of a and along a conjugate path, in contradiction to the general uniqueness argument above. 
(Note that solutions at a = i and a = —i should be considered equivalent; the change in 
sign of OL simply means that the path is traversed in the opposite direction.) Now it is 
clear that there can be no integral curve from a; = to a; = A, because it would have to 
pass either over or under x = 1, and thus could not lie entirely along the real axis. Thus, 
there must be exactly two integral curves (since this is the minimum possible), from to 
1 (which occurs at a = 1) and from 1 to A (at a = i). These correspond to real and 
imaginary 3-spheres in the Calabi-Yau manifold, intersecting transversely at a point. 
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Next, consider interpolating between these configurations by holding two of the roots 
fixed and moving the third root along a line midway between them. We do this by varying 
A in Eq.( |5.3| ) from ^ + ^ 2 ^° I' ^lo'^S ^^i^ ^^^^ Re(a;) = |. We find numerically that at 
approximately A = ^ + .231i, the number of integral curves jumps from 3 to 2. The curve 
of marginal stability in the A-plane has roughly the appearance shown in Figure 7, and is 
preserved by the modular transformations 

A — ^i A — ^1-A, A — ^A. (5.4) 

A 

In each of the regions containing the real A-axis, the number of stable BPS states is 2; 

elsewhere in the A-plane, there are 3 such states. By a local analysis, it can be shown 

that the components of the curve intersect the real axis at angles of ±60°, and approach 

asymptotes also making angles of ±60°. 




Fig. 7; Curves of marginal stability for n = 3. 

The general An singularity with P{x) an nth order polynomial may be analyzed by 
the same methods. Again, when the roots are arranged with Z„ symmetry, we will show 
that the maximum possible number n(n — l)/2 of stable BPS states is realized. And when 
the roots are colinear, the number of stable BPS states is minimized. Indeed, in the latter 
case, the same argument we gave for n = 3 extends to show that there are precisely n — 1 
integral curves, connecting the roots in sequence, with a alternating between 1 and i. 

On the other hand, suppose that the roots are located at the nth roots of unity, 
corresponding to the polynomial 

P{x) = x'^ - 1 , 
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We will now show that there is a unique integral curve connecting any two roots, for a 
total of n{n — l)/2 stable hypermultiplet states. 

Assume for simplicity that n is even. In general, all integral curves of Eq.( p.l| ) for an 
nth-order P{x) that extend out to infinity in the x-plane must asymptote to one of n + 2 
lines of constant phase. This is because, for large x, Eq. (|5.1| ) has the asymptotic solution 



x{t) = {{n + 2)at) 



l/(n+2) 



As we have discussed, there are exactly 3 integral curves emanating from each root; each 
of these must end either at another root or at one of these n + 2 asymptotic infinities. We 
consider the graph formed by the set of all integral curves starting or ending at a root, 
such as the graph shown in Figure 8 for the case n = 6 and a = 1. 




Fig. 8; Graph of integral curves ending on roots for n = 6 and a = 1. 



An important fact, which severely constrains the topologies of allowed graphs of this 
type, is that no two integral curves from a single root may approach the same asymptotic 
infinity. The proof is almost identical to the above argument that no two roots can be joined 
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by two distinct integral curves. If this did happen, we could consider the concatenation of 
two integral curves between a given root and a given infinity. A small deformation would 
produce a curve beginning and ending at the same infinity, with net index of +1/2. But 
as we have seen, closed curves always lead to indices of or less, a contradiction. 

We now argue that, for a=l, there are integral curves connecting complex-conjugate 
nth roots of 1, and all other curves in the graph are unbounded. To see that only conjugate 
points can be connected, suppose that two non-conjugate roots Xj = exp(27rzj/(n + 2)) and 
Xj+k = exp(27rz(j+/c)/(n+2)), were joined by an integral curve. Then complex conjugation 
would map this curve into a different integral curve, also with a=l, between a pair of 
roots conjugate to the original pair. By Z^ symmetry, at other values a = exp{2nik/n), 
we obtain n rotated pairs of integral curves. Generically, this procedure leads to multiple 
geodesies between Xj and Xj+fc, at distinct values of a. (A special case, when the original 
two roots are related by refiection about the imaginary axis, can be ruled out by showing, 
as we do below, that integral curves joining such points already exist for a = i.) 

Next, we need to show that, for a = 1, all conjugate roots are indeed connected 
by integral curves. We first assume that n is even. The roots x = 1 and x = —1 are 
special; they are not connected to any other roots, and each of them has 3 integral curves 
which must terminate at 3 different asymptotic infinities, as seen in Figure 8. By complex 
conjugation symmetry, one of the integral curves from x = 1 lies entirely along the real 
axis; the other two are asymptotically parallel to the lines x{t) = exp(±27rz/(n + 2))i. In 
traversing a very large circle clockwise in the x-plane, the winding index of any vector 
field defined by (|5.1|) will be (n + 2)/2. As we traverse the circle clockwise from the point 
where the first integral curve from x = 1 crosses this circle to the point where the third 
curve crosses, the index shifts by at least +1. Now consider any root connected to its 
conjugate; it must also have two unbounded integral curves, which similarly contribute at 
least +1/2 winding to the total index. If all conjugate roots are connected, the total index 
will be at least {n + 2)/2; if any are not connected, their individual contributions to the 
index will be at least +1, and the lower bound on the index will be greater than its actual 
value. Therefore, for a = 1, all conjugate roots are connected, and we have thus shown 
the existence of ^ — 1 integral curves of finite length. The graph of integral curves for the 
case n = 8 is shown in Figure 8. 

The same sort of argument holds for each a which is an nth root of 1 . In each of the 
n/2 cases {i.e., not distinguishing between a and —a), we can exploit the symmetry about 
the line at to obtain f — 1 integral curves, for a total of ^(^ — 1) distinct cases. We can 
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also make use of the symmetry about the hue at when a is a (2n)th root of 1; this gives 
an additional (§)^ integral curves, for a total of n{n — l)/2. Finally, if n is odd, a similar 
counting scheme works: each nth root of 1 corresponds to a distinct symmetry axis, and 
each such value of a leads to (n — l)/2 integral curves, for a grand total of n{n — l)/2. 
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